A new dissipation term for finite-difference simulations in Relativity

We present a new numerical dissipation algorithm, which can be efficiently used in combination with centered finite-difference methods. We start from a formulation of centered finite-volume methods for Numerical Relativity, in which third-order space accuracy can be obtained by employing just piecewise-linear reconstruction. We obtain a simplified version of the algorithm, which can be viewed as a centered finite-difference method plus some 'adaptive dissipation'. The performance of this algorithm is confirmed by numerical results obtained from 3D black hole simulations.Comment: Talk presented at the Spanish Relativity Meeting (Tenerife 2007

New Formalism for Numerical Relativity

We present a new formulation of the Einstein equations that casts them in an explicitly first order, flux-conservative, hyperbolic form. We show that this now can be done for a wide class of time slicing conditions, including maximal slicing, making it potentially very useful for numerical relativity. This development permits the application to the Einstein equations of advanced numerical methods developed to solve the fluid dynamic equations, {\em without} overly restricting the time slicing, for the first time. The full set of characteristic fields and speeds is explicitly given.Comment: uucompresed PS file. 4 pages including 1 figure. Revised version adds a figure showing a comparison between the standard ADM approach and the new formulation. Also available at http://jean-luc.ncsa.uiuc.edu/Papers/ Appeared in Physical Review Letters 75, 600 (1995

Gowdy waves as a test-bed for constraint-preserving boundary conditions

Gowdy waves, one of the standard 'apples with apples' tests, is proposed as a test-bed for constraint-preserving boundary conditions in the non-linear regime. As an illustration, energy-constraint preservation is separately tested in the Z4 framework. Both algebraic conditions, derived from energy estimates, and derivative conditions, deduced from the constraint-propagation system, are considered. The numerical errors at the boundary are of the same order than those at the interior points.Comment: 5 pages, 1 figure. Contribution to the Spanish Relativity Meeting 200

CP Asymmetries in B to f_0 K_S Decays

We consider the branching ratio and the CP asymmetries in B to f_0(980)K_S decay to the end of determining the deviation of the time-dependent CP asymmetry from sin(2 beta) arising from Standard Model physics. We obtain Delta S_{f_0 K_S} within the context of the QCD factorization framework for the B to f_0(980)K_S decay amplitudes assuming the f_0(980) is a q\bar{q} state and employing a random scan over the theoretical parameter space to assess the possible range in Delta S_{f_0 K_S}. Imposing the value of the experimental branching ratio within 1 sigma and 3 sigma, respectively, of its central value as a constraint, we find the range of Delta S_{f_0 K_S} to be [0.018, 0.033] for a scan in which the parameters are allowed to vary within 1 sigma of their central values and the range [-0.019, 0.064] for a scan in which the parameters vary within 3 sigma of their central values.Comment: 27 pages, 10 figures, references adde

Observation of B+ ---> a(1)+(1260) K0 and B0 ---> a(1)-(1260) K+

We present branching fraction measurements of the decays B^{+} -> a1(1260)^{+} K^{0} and B^{0} to a1(1260)^{-} K^{+} with a1(1260)^{+} -> pi^{-} pi^{+} pi^{+}. The data sample corresponds to 383 million B B-bar pairs produced in e^{+}e^{-} annihilation through the Y(4S) resonance. We measure the products of the branching fractions: B(B^{+}-> a1(1260)^{+} K^{0})B(a1(1260)^{+} -> pi^{-} pi^{+} pi^{+}) = (17.4 +/- 2.5 +/- 2.2) 10^{-6} B(B^{0}-> a1(1260)^{-} K^{+})B(a1(1260)^{-} -> pi^{+} pi^{-} pi^{-}) = (8.2 +/- 1.5 +/- 1.2) 10^{-6}. We also measure the charge asymmetries A_{ch}(B^{+} -> a1(1260)^{+} K^{0})= 0.12 +/- 0.11 +/- 0.02 and A_{ch}(B^{0} -> a1(1260)^{-} K^{+})= -0.16 +/- 0.12 +/- 0.01. The first uncertainty quoted is statistical and the second is systematic

Observation of the semileptonic decays B ---> D* tau- anti-nu(tau) and evidence for B ---> D tau- anti-nu(tau

We present measurements of the semileptonic decays B- --> D0 tau- nubar, B- --> D*0 tau- nubar, B0bar --> D+ tau- nubar, and B0bar --> D*+ tau- nubar, which are potentially sensitive to non--Standard Model amplitudes. The data sample comprises 232x10^6 Upsilon(4S) --> BBbar decays collected with the BaBar detector. From a combined fit to B- and B0bar channels, we obtain the branching fractions B(B --> D tau- nubar) = (0.86 +/- 0.24 +/- 0.11 +/- 0.06)% and B(B --> D* tau- nubar) = (1.62 +/- 0.31 +/- 0.10 +/- 0.05)% (normalized for the B0bar), where the uncertainties are statistical, systematic, and normalization-mode-related

Search for CP violation in the decays D0 ---> K- K+ and D0 ---> pi- pi+

We measure CP-violating asymmetries of neutral charmed mesons in the modes D0 --> K- K+ and D0 --> pi- pi+ with the highest precision to date by using D0 --> K- pi+ decays to correct detector asymmetries. An analysis of 385.8 fb-1 of data collected with the BaBar detector yields values of aCP(KK) = (0.00 +/- 0.34 (stat.) +/- 0.13 (syst.))% and aCP(pipi) = (-0.24 +/- 0.52 (stat.) +/- 0.22 (syst.))%, which agree with Standard Model prediction

Improved estimates of rare K decay matrix-elements from Kl3 decays

The estimation of rare K decay matrix-elements from Kl3 experimental data is extended beyond LO in Chiral Perturbation Theory. Isospin-breaking effects at NLO (and partially NNLO) in the ChPT expansion, as well as QED radiative corrections are now accounted for. The analysis relies mainly on the cleanness of two specific ratios of form-factors, for which the theoretical control is excellent. As a result, the uncertainties on the K+ --> pi+ nu nubar and KL --> pi0 nu nubar matrix-elements are reduced by a factor of about 7 and 4, respectively, and similarly for the direct CP-violating contribution to KL --> pi0 l+ l-. They could be reduced even further with better experimental data for the Kl3 slopes and the K+l3 branching ratios. As a result, the non-parametric errors for B(K --> pi nu nubar) and for the direct CP-violating contributions to B(KL --> pi0 l+ l-) are now completely dominated by those on the short-distance physics.Comment: 16 pages, 1 figure. Numerical analysis updated to include the recent Kl3 data. To appear in Phys. Rev.

Formulations of the 3+1 evolution equations in curvilinear coordinates

Following Brown, in this paper we give an overview of how to modify standard hyperbolic formulations of the 3+1 evolution equations of General Relativity in such a way that all auxiliary quantities are true tensors, thus allowing for these formulations to be used with curvilinear sets of coordinates such as spherical or cylindrical coordinates. After considering the general case for both the Nagy-Ortiz-Reula (NOR) and the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formulations, we specialize to the case of spherical symmetry and also discuss the issue of regularity at the origin. Finally, we show some numerical examples of the modified BSSN formulation at work in spherical symmetry.Comment: 19 pages, 12 figure

Measurements of Partial Branching Fractions for anti-B ---> X(u) l anti-nu and Determination of |V(ub)|

We present partial branching fractions for inclusive charmless semileptonic B decays Bbar --> X_u ell nubar, and the determination of the CKM matrix element |V_{ub}|. The analysis is based on a sample of 383 million Y(4S) decays into B Bbar pairs collected with the BaBar detector at the PEP-II e+ e- storage rings. We select events using either the invariant mass M_X of the hadronic system, the invariant mass squared, q^2, of the lepton and neutrino pair, the kinematic variable P_+ or one of their combinations. We then determine partial branching fractions in limited regions of phase space: Delta B = (1.18 +- 0.09_{stat.} +- 0.07_{sys.} +- 0.01_{theo.}) x 10^{-3} (M_X 8 GeV^2/c^4). Corresponding values of |V_{ub}| are extracted using several theoretical calculations